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Abstract 



We present new experimental results on the dielectric constant in orientationnal glasses 
Ki_^ Li^Ta03 (KLT) with x < 0.05, together with the detailed (analytic and numerical) study 

of a model which attributes the observed aging to the motion of the walls of polarization 
domains. We show that the dielectric constant after a positive temperature jump goes through a 
maximum as a function of the subsequent time. This observation and those previously reported 
(aging, cooling rate dependences, etc) are compared with the predictions of the model, in which 
the variations of the dielectric constant are attributed to the change of polarization domain wall 
area. The total area decreases by domain growth and increases by nucleation of new small 
domains inside the large ones. These two opposite variations are both hindered by static 
random fields (equivalent to energy barriers) due to the frozen dipoles borne by the off-center 
Li+ ions. Many results are well explained by a model with a single energy barrier. However, 
some effects can only be understood if a broad distribution of energy barriers is assumed. We 
use the experimental data to determine this distribution and find it to be unimodal with a width 
comparable to its most probable value. 
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1. Introduction 



Ergodicity is a usual assumption of Statistical Mechanics. A system that satisfies this 
hypothesis is able to visit all allowable points of its phase space which are specified by the 
macroscopic constraints. In other words, the time average and the ensemble average of any 
thermodynamic quantity coincide. This general rule encounters some important exceptions: for 
instance, when a phase transition occurs, the associated symmetry breaking induces ergodicity 
breaking into mutually inaccessible regions of phase space [1]. 

Indeed, the dynamics and time scales play an important role. If the phase space is 
entirely visited in a finite time, the system is ergodic. If the phase space is splitted into several 
sub-spaces separated by barriers which diverge with the size of the system, the representative 
point of the instantaneous state may not be able to visit all space, even after times which 
diverge, thought not fast enough to allow it to surmount the diverging barriers. Then it is said 
that ergodicity is really broken (true ergodicity breaking). Another example is a case in which 
the arrival point (free energy minimum) in the phase space is separated of the departure point by 
an infinite number of finite barriers. Then the representative point will certainly arrive into the 
minimum, but this may take an infinite time (weak ergodicity breaking [2]). 

From the experimentalist's point of view, things are not so clearly cut: actually, it may 
be practically impossible to distinguish between a finite time of one year and an infinite time. 
This somewhat fuzzy case (effective ergodicity breaking) corresponds to slow dynamics as 
observed in spin-glasses and many other systems. 

The slow dynamics of orientational glasses (OG) has been extensively studied in the last 
years. The attention has been focused on two points: aging and ergodicity breaking. OG are 
materials which bear electric or elastic moments with orientational degrees of freedom 
distributed on random sites [3, 4]. In these crystals, each dipole or quadrupole is submitted to 
conflicting interactions from the Ugand field and the fields generated by other moments located 
at random. Since the bearers of the moments are sticked to a given site, the disorder of their 
interactions is frozen in. This leads to frustration: the dipoles are not free to rotate; their 
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moment, locked by high potential barriers, are not able to minimize the total free energy of the 
system. 

Among other OGs, the potassium-lithium tantalate series K2_j^ Lij^Ta03 (KLT in 
short) has been particularly studied. The time-dependent response after different perturbations 
was measured. Examples of such response are: the birefringence after switch off of a 
mechanical stress [5], the elastic constant c(co,t) [6] and, mainly, the dielectric constant e(co,t) 
after different thermal histories (regular temperature ramp, quenching, temperature jump, 
temperature cycle, ...) where 0) is the circular frequency related to the frequency / by 
(0 = 2% f [7-9]. The alternative dielectric constant £(co,t) was found to slowly decrease with 
time, showing a non stationary behavior or aging. Using a fitting function, an extrapolated 
"asymptotic" value £(0),oo) has been deduced. The dependence of e(G),oo) on the cooling rate 
was discussed in terms of effective ergodicity breaking. In the end, a tremendous amount of 
data was thus gathered. But no microscopic model was proposed in order to explain the 
measurements. 

In this paper we present a model which is intended to describe the whole set of data 
recorded on a series of KLT crystals. The model is based on the idea that the time dependent 
part of the dielectric constant (and to any other susceptibility such as the elastic constant) is due 
to the contribution of the wall of polarization domains, the size of which slowly evolves 
towards their equilibrium value. A short report has been recently published [10]. Concerning 
our experimental data, many of them are already published [6-9] and are only briefly recalled. 
Others are new and have a strong heuristic role for improving the model; therefore, they are 
reported in more detail. 

The theoretical part of the paper is organized along two main lines: i) we first describe 
the scaling arguments, based on a domain growth picture, which lead to explicit expressions for 
the evolution of the dielectric constant for typical experiments; ii) we then perform a few 
numerical computations, on some particularly demonstrative experiments, which check or 
extend the validity of the previous calculations. In both cases, our attention is focused on the 
long time behavior. 
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We must emphasize that our model is fundamentally ergodic. Nevertheless, because of 
the time scales in play, it presents many of the features of ergodicity breaking: this is effective 
ergodicity breaking. 
2. Experiments 

2.1. The Material - Pure KTa03 (KT) is a cubic perovskite crystal. It is called an 
incipient ferroelectric because its trend towards ferroelectricity is prevented by quantum 

fluctuations. Upon cooling, the displacements of the ions slow down but never become 
ordered. However, the correlation length of the polarization steadily increases as shown by the 
dielectric constant e which strongly rises as the temperature tends towards K. The symmetry 

of the soft modes responsible for the enhancement of £ corresponds to the displacements of the 

Ta^+ ions along one of the six [100] directions [11, 12]. 

On random substitution of Li+ ions to K+ ions in a KTa03 crystal, the tendency 

towards ordering of the Ta5+ ions is reduced [13]. Because of their small size, the Li+ ions take 
off-center positions, referring to the centro-symmetric site. This generates electric dipoles at 
random sites which are aligned along one of the six [100] directions. Every dipole interacts with 
the others through the polarizable surrounding medium. At high temperatures, the dipoles 
constantly reorient between their six possible positions and the crystal is in a paraelectric phase. 
As temperature is lowered, the dipole reorientations slow down: the characteristic time for a 
jump over the barrier of Tg = 1000 K is longer than 3 years for T < 24 K. According to the 
Li+ concentration x, the dipoles gradually freeze into a disordered phase (x < 0. 022) or they 
suddenly freeze into a cooperative order (x > 0.022). Indeed, for any concentration x larger 
than 0.01 there is a transition temperature which separates the paraelectric phase from either 

the glassy or the ferroelectric phase. The phase obtained for low concentration is a frozen 
configuration with neither long-range ferroelectric order nor static polarization; it is an 
orientational or dipolar glass. We emphasize that the phase diagram thus described concerns the 
state of the Li+ ions, not the state of the tantalate lattice that is always in a paraelectric phase. At 
low temperatures, its dielectric constant 8, while smaller than in pure KT, is rather large; this 
indicates the presence of large domains of fluctuating polarization. 
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2.2. The Samples - We have studied a series of mixed single crystals of formula 
^l-x Lijt:Ta03. In addition to a nominally pure sample (x < 0.001), we have examined five 
samples, the concentration x and the transition temperature Tf^ of which are listed below: 0.0 11 

(27 K); 0.017 (35 K); 0.025 (46 K); 0.030 (52 K); 0.050 (73 K). The lithium concentration 
X is deduced according to the empirical law r^^(K) = 535 x ' proposed by van der Klink et 

al. [13]. The error bars are Ax = ± 0. 001 . 

2.3. The Measurement Methods - All the faces of the samples are oriented perpendicular 
to a [100] direction of the cubic crystal. With chromium electrodes on two opposite end-faces, 
they form capacitors. The real and imaginary parts of the dielectric constant 
e{(0) = e\(0)-i e"{(0) are related to the capacitance and loss. They were measured with a 
Hewlett-Packard 4 192 A impedance analyser at seven frequencies ranging from 1 kHz to 
1 MHz. Some elastic constant measurements were also performed, in the 100 MHz range. 
They provide the real and imaginary parts of the elastic constant c{(0) = c'(co) - i c"(co). In 
KLT the imaginary parts are much smaller than the real part and they are not used in what 
follows. It was checked that the response is linear. All the measurements were performed 
without any applied static electric field or mechanical stress, i.e., without any symmetry 
breaking. Since the capacitance C and the real part e' of the dielectric constant are simply 
proportional through a geometrical factor, our data are generally displayed in pF in the figures 
as they were recorded. 

2.4. Results Already Known — Our experiments can be classified in several types 
according to the procedure used to attain the experiment temperature T^^p. 

(i) Variable cooling ramp. - The sample is regularly cooled from a high temperature through the 
transition temperature T^^ down to T^^p at a constant rate R = dT/dt . Then, the dielectric 

constant e'(co, t) starts from e\co, 0) and steadily decreases. This aging is very well fitted by a 
power law function £'(©,?) = e(©,°°) + Ae + ?o)Ao) " for times as long as 400000 s. 
However, the main result of this analysis is that the apparent "asymptotic" value £{(0,°°) 
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depends on the cooling rate with a larger | R | leading to a higher £(co, o°). More precisely, the 
dependence is approximately e(co,oo) = p + q log| R | (see figures 5 and 6 in reference [8]). All 
the samples share this behavior with a positive q, except for the nominally pure sample for 
which q is slightly negative. 

(ii) Interrupted cooling ramp. - The sample is cooled at a constant rate R down to the 
quenching temperature Tq below Tf^., from where it is very rapidly cooled at the rapid rate 

Rg =15 R down to Tg^p. The same R and Rg are used for a series of quenching temperatures 

with Tg^p <Tq < Tff.. Figure 11 in reference [8] shows the dependence of the initial value 

e((O,0), at the beginning of the isothermal evolution at Tg-^p together with the "asymptotic" 

value e(co, o°) for the sample with x = 0.025. They both have the same behavior: they depend 
on Tq if it is close to , but they become independent of Tq if Tq is small enough. The 

comparison of our samples with x = 0.025 and x = 0.011 shows that their behaviors are 
correlated with the dependences on the cooling rates: the larger coefficient q ma cooling ramp 
experiment corresponds to the larger curvature in a quenching experiment. Therefore, these two 
results demonstrate that what is important for the further evolution of e'((0, t) is the manner in 

which the transition temperature is crossed. 

(iii) Temperature cycles. - Two temperatures Ti and T2 are used in three steps. After the 
sample was regularly cooled down to Ty, it is kept in isothermal conditions during ti. Then, 
the temperature is suddenly changed to T2 and kept constant for ?2- Finally, the sample is 
suddenly driven back to T^. Measurements are continuously taken all along the three periods. 
The temperature T2 may be smaller or larger than the temperature Ty. If T2 > 7i, the sojourn at 
T2 acts as a simple annealling [8]. If T2 <Ti, the experiment is much more interesting since a 
paradox appears: on the one hand, the system ages at T2 during the time t2 ; on the other hand, 
after the return to Tj, a continuous curve joining the first and third periods is obtained by 
simply translating the third period, just as if the second period had not existed (see figure 14 in 
reference [8]). Remarkable transient effects were observed after positive temperature jumps. 
They are reported in detail in the next sub-section. 

Similar experiments have been performed in spin-glasses (for a review, see [14]). 
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2.5. New Results - During the course of temperature cycle experiments we had observed a 
surprising behavior which immediately followed a positive temperature jump. At the jump the 
dielectric constant starts from a high value at the low temperature Ti and rapidly decreases to a 
smaller value at T2>Ti. But e' first reaches a low value, then it immediately increases back, 
passes by a maximum and finally resumes its slow downward aging [9]. We present here a 
more complete report on this phenomenon studied in the sample with x = 0. 017. 

At fixed low temperature Ti and for a fixed duration ti, we have performed jumps to 
different high temperature T2>Ti. The larger the temperature jump Tj-Ti, the larger the 
downward jump of e'. The behavior of e'(co,t) is always non monotonic (see figure 1), 
whatever T2 , lower or higher than the transition temperature Tff. . In this last case the final 

decrease is very small and rapid. 

At fixed temperatures Ti and T2 , we have studied the role of the duration ti which was 
varied between 500 s and 80000 s. The curve corresponding to the shortest ti is the steepest 
and it encompasses all the other. For the longest ti, the curve has a flat maximum, lower than 
all the other maxima. When ti increases, the location of the maximum is shifted towards larger 
times. For long times, all curves merge into a common curve (see figure 2). We have analysed 
the curves as follows. They are decomposed in two parts: a long time decrease, shared by all 
the curves and well described by a single power law and a short time increase. The difference 

between the experimental data and the best power law separates the initial increase that may be 
described as an exponential return 4ej„^^ [l - exp(-r/Tj„(,^)]. The amplitude 4ej„^^ rises with 

ti but saturates for large ti. Surprisingly, ^ej>jcr is smaller for larger T2 (see figure 3). The 
characteristic time T,^^^ is nearly independent of ti (see figure 4) and it is longer for lower 
temperature T2 , approximately following the Arrhenius law with an activation barrier energy 
T}j of the order of 15 K. Its value Tj^^^ = 1200 s is short in comparison to that of aging in 
long experiments (of the order of ti). 

2.6. Dielectric Measurements and Elastic Measurements - The results obtained by 
dielectric measurements and elastic measurements show unexpected similarities. After a regular 
temperature ramp, the real part of the dielectric constant e'(co, t) and the real part c\g), t) of the 
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elastic constant both isothermally decrease as time flows. After a positive temperature jump, 
both £\(0,t) and c\0),t) first increase, then pass by a maximum and finally decrease [6, 8]. 
This common behavior is surprising since it is well known that the interaction of a wave with 
relaxing entities always slows down the wave. This is true for the electromagnetic wave 
velocity v^^ as well as for the elastic wave velocity Vg/. Therefore, since = q l^£' fi and 
^el ~ '\I^'/P ( is lig'^t velocity in vacuum, fj. is the magnetic permittivity and p is the 
specific mass) the dielectric constant e' and the elastic constant c' are expected to present 
opposite behaviors. This is actually observed in KLT in the temperature range 40-60 K where 
e' increases while c' decreases, due to the relaxation interaction with the Li+ ions jumping 
over their potential barriers. Therefore, the two quantities which must be compared are the 
dielectric constant e' and the elastic compliance s' = l/c' , inverse of the elastic constant or 
stiffness (the tensorial character of these quantities is neglected here). In this context, the similar 
behaviors observed on e' and c' during aging were not expected at all. 

The preceding remarks rule out, as source of aging, any relaxation process such as the 
thermally activated jumps of the Li+ ions between equivalent off-center positions. They prompt 
to look for new interaction processes such that e' and c' must behave in the same manner. 

3. The Model 

We recall here some important features concerning KLT which are relevant for the 
model. Firstly, the dielectric constant is large and increases as temperature is lowered, thus 
indicating the presence of electric polarization domains of the host lattice, oriented along one of 
its six [100] directions, with growing size. Secondly, below 40 K the Li+ ions are completely 
frozen in and the dipoles they bear, also oriented along one of the six [100] directions of the 
lattice, generate electric static random fields. Thirdly, the highly polarizable host crystal is 
locally sensitive to the presence of the Li+ ions: around any of them there is a distortion of the 
lattice which gives it a large effective dipolar moment, as measured in [15] and explained in 
[16]. 
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Figure 5 shows the results of a numerical simulation of a three-dimensional assembly 
of Ising spins with ferromagnetic interactions in a random magnetic field. We present a two- 
dimensional cross-section of the system at four different times. The up-spins are represented by 
black dots and the down-spins by white dots. These pictures provide some information on the 
geometry of the domains of ordered spins. Indeed, even if their shapes are very whimsical, 
some features are clearly seen: they have a typical size which increases with time and there is 
also a visible size dispersion. 

These facts are in favor of the role of the polarization domains in the time-dependent 
behavior of the dielectric constant e' . Since e'(t) decreases when the typical domain size r(t) 
increases (because the volume fraction of the walls then decreases), we will assume in the 
following that the major contribution to the aging part of the dielectric constant comes from the 
domain walls, and not from the bulk. 

Actually, the pictures of figure 5, shown as an illustrative example for the domain 
shape, may be misleading and we have to emphasize the physical differences between those 
Ising spin domains and the polarization domains in KLT. On the one hand, the Ising spin 
domains which bear a quasi static magnetization are dipolar domains. On the other hand, the 
KLT lattice is in a paraelectric phase of an incipient ferroelectric. This means that there are polar 
soft modes, the frequency of which lowers when the temperature decreases, as evidenced by 
the increase of e' . These domains are the forerunners of the aborted ferroelectric transition at 
K. At any time they bear an electric polarization which oscillates at the soft mode frequency. 
Hence, the mean polarization taken over time intervals longer than 10~^^ s (as for our 
experiments), vanishes. Consequently, for our purpose the KLT domains are rather 
quadrupolar than dipolar, the two directions along any of the three cubic axes cannot be 
distinguished and, finally, there are only three types of domains in KLT, oriented along the 
three [100] axes. An orientable entity in KLT passes from a domain to a neighboring one by a 
71/ 2 -flip while for Ising domains this needs a 7t-flip. 

3.1. First Insight on the Model - We assume that the time-dependent contribution to the 

complex dielectric constant £ (and to the complex elastic constant c too) is due to the 
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polarization domain walls. These walls change in two manners. They have an oscillating 
motion, forced by the applied alternative electric field (or elastic strain); the motion is not in 
phase with the field and the subsequent hysteresis induces loss and dispersion. Moreover, if the 
domains are out of equilibrium, the walls also have a natural slow motion towards the 
equilibrium size, defined as the correlation length ^(T) of the fluctuations of polarization. The 
length (§(r) depends on temperature; its square is proportional to the equilibrium value of the 
dielectric constant, and consequently, ^(T) = ^f^-^jY^^Jf where the correlation length 

^ff. = ^{Tfy.) at the transition temperature Tjj. is taken as reference. The natural motion of the 

walls is hindered by the static random fields generated by the Li+ ions when they are frozen 
(for T < Tff.); the motion is free when the Li+ ions are permanently re-orientating (for T > Tf^), 

producing zero field in average. As a result of the two movements, the complex susceptibilities 
measured at the circular frequency (0 also depend on time and read e(co,t) and c(co,t). We 
notice that the oscillating motion of the domain walls corresponds to polarizability and elasticity 
added in parallel to those of the lattice, as discussed in sub-section 2.6. Because most 
experiments were dielectric measurements, in what follows we focus our attention on e(co,t) 
only. 

Moreover, we assume that the loss and dispersion are proportional to the total wall area. 
The proportionality factor is an efficiency coefficient which contains the cosine of the angle 
between the field and the normal direction of the wall at any point of the surface. 

As a first step, we assume that the size dispersion is weak enough to be neglected. Then 
all the domains have the same size r(t). Then the area of the surface of a domain is <7(t) ~ r(t) 
(a numeric coefficient of geometric origin is omitted), while the number of domains in the 
sample of volume Vis n(t) ~ v/ r(t)^ (here too, a numeric coefficient is omitted). Finally, it is 
found that the total wall area is E{t) = C V/r(t), where the coefficient C is a dimensionless 
geometrical factor which takes into account the orientation of the field and a factor of 1/2 
which prevents to count twice the surfaces which are shared by contiguous domains. In most 
experiments, isothermal aging is studied after cooling. Then the size r(t) at temperature T is 
smaller than ^(T) = ^f^-^jT^^Jf and the trend of the domains towards equilibrium can only be 

growth. If r(t) increases, then Z(t) and £'(Q},t) decrease, as actually observed. 
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Pinning impurities strongly impede the domain wall motion and lead to thermally 
activated motion [17]. We model this by writing a simple rate equation for the isothermal 
evolution of their size p(t) : 

^^^^ (1) 

dt T(p,r) 

where 

t(p,T)=t^ exp 



T 

6 is an exponent describing how the energy barriers typically grow with the length scale p 
involved in the motion, T(j measures the typical strength of the pinning field while a is of the 
order of the lattice cell parameter. The exponent 6 depends on the coupling between the polar 
entities: 9 = 4/3 for nearest neighbor interactions while 6 = 2 for long range dipolar 
interactions. In the following we use 6 = 1, a value not far from the previous realistic ones, 
which allows equation (1) to be solved exactly as: 

p(t) = Po+a(T/T^)\n[l + t/T(po,T)] , 

where pQ is the initial value of the domain size. This function increases rapidly as 
p(t)~pQ + a(T/Tfj)t/T(pQ,T) for t/T(pQ,T)«l and slowly as 
p(t) ~ a (T/T(y) ln[?/T(Po,r)] for t/T(pQ,T) » 1. Accordingly, the free size p(t) can go to 

infinity. In fact, the true size r(t) is necessarily limited by the thermal correlation length ^(T), 
then 

r(t) = mm{p(t),^(T)} , 
and the corresponding relaxation time becomes 

fr-T ^ 

T(rj-,r)=T^exp ^-f with ri = r(ti) . 
{ aT J 

The relaxation time has the Arrhenius form with the barrier height T^j, measured in Kelvins; 
Too is a microscopic time. This form means that the activation process responsible for the 
evolution for T <Tfj. is limited by the effective barrier height T(j(ri) = (ri/a)T(^ which 
depends on the initial domain size rj scaled by the lattice cell parameter a. If T > the Li+ 
ions are not frozen, then Tq- = and T(rj,T) = z^o is very short. 
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The simplest experiment begins with a controlled cooling at constant cooling rate 
R = dT/dt from the initial temperature r,„,Y above Tf^. to the experiment temperature T^^p 

below Tff.. The temperature variation reads T(t) = 7^^ + Rt with R<0. The relevant cooling 
time is =(Tf^- Tg^p)l\R\. The correlation length depends on time through T(t). It reads 
<^(0 = ^fr^jTir/(Tfr + Rt) for 0<t<tc (or T^^p <T< Tf^) if the time origin is taken at the 
instant when Tf^. is crossed. The evolution velocity of ^(t) and the domain growth velocity 
respectively read 

dt 2 Tit) d< T(p,T{t)) 

As long as p > the domains are able to reach their maximum size ^(T) at any time during 
cooling. This is no longer true for t>t* where the time t* is determined by the condition 
p(t ) = ^(T(t )). Two cases have to be considered. 

If t* >t(., during all the cooling and afterwards, the domains are in equilibrium. The 

domain size is 

r{t) = ^{T{t)) for < r < , 

r{t)=^{T^^) for t>tc . 

Consequently, there is no isothermal aging since r(t) does not vary for t>t^. 

If <t* <t^, the domains are in equilibrium up to t* and fall out of equilibrium when 
t>t . Afterwards, they non isothermally grow from r =T{t ), where r{t ) = q\Tfj. + Rt I, 

down to Tg-j^p and finally, isothermally grow. Therefore, the domain size is 

r{t) = ^{T{t)) fox Q<t<t* , 



for t* <t <tc , ( 3 ) 



c ' 



J df 

r{t) = min [r^ + a (T^xp/Ta) ln[l + (t- ?c)/t(^c. Texp)] , ^(Texp)} for t > 

with = r{tr). 



The values of R and T^y determine the case to be considered. Indeed, the time t is the 
solution of the equation 
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dt 



which reads, with T =Tff. + Rt , 



2 T* 



- a 



(T*/T^) 



r(^(r*),r*) ■ 

In this equation the right-hand side is the derivative of p(t) written with t = t* and 
r(t*) = ^(T*) as initial values. After a trite calculation we find 

x2/3 



r = 



Jr. 

a 



tr 



1/2 



In 



1^1 ^oo ^tr T, 



1/2 



tr 



(4) 



The logarithm is a slowly varying function if its argument is much larger than 1. Therefore the 

unknown temperature T* may be replaced in the right-hand side by any temperature between 
Tf^ and Tg^p without noticeable change. We arbitrarily choose T* = T^^ and set 



2 a r,,2 



T^oo ^tr 



to find, as a good approximation. 



/ e 7^ 1/2 7^ 



a ln(i?b/|^|) 



With z-^=10~^^s, i? = -0.1K.s~^ and r^^ = 35K we obtain r*=4.7K if 
T(j <^f^/(3 = 50 K or r* = 37 K if T^j ^trl^ ~ 1^00 K, while an order of magnitude for Rq is 
Rq ~ 10^^ K.s~\ which is extremely large. Therefore, in the following, we shall often assume 
that lni?o » In the first case T* <Tg^p = 5K and the domains are always in 

equilibrium: the characteristic time T is still of the order of the microsecond at the lowest 
temperature Tg^p = 5 K. In the second case T* > T^j., the domains fall out of equilibrium as 

soon as T^^ is crossed (above T^^ they have the equilibrium size (^(7)): the characteristic time 
is already T = 30000 s at T^j. and longer than the age of the universe at Tg.^p. In both cases 
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there is no aging at T^^p, either because the domains are too rapidly in equilibrium or because 

they are completely frozen. Therefore, aging must be looked for in the intermediate range. 

If T* is close to the difference between the coherence length and the domain size 

can be expanded in powers of {T* - T^^) and to second order it reads 

^{t^)-p{tc) = ^tr ^''ZK .5/3 On(^/l^|)f [5 + 3 MRo/\R\)] (5 ) 





* 


2 




T -T 

^ ^ exp 





8 r,//3 7^^5/3 



htr 
V a J 



Writing hi(i?o/|^|) = In(^) [l- ln(|i?|)/ln(i?o)] and using ]n(\R\)/\n(RQ) «l, itis clear that 
the result may indeed be put in the form ^(Tg^p) - p(t^) = A + B ln(|i?|)/ln(i?o)- 

The corresponding variation of e' is 

e'(rJ-e'(oo) = C + Z)ln(|i?|)/ln(i?o) • (6) 

as experimentally observed [8]. 

Another simple experiment is quenching. It begins at the temperature Tf^, with a 
controlled cooling at constant rate R = dT/dt , interrupted at the temperature by a sudden 
(instantaneous, in principle) jump down to T^^^p. The quenching happens at 
tq = (Ttr - Tq)/\R\. The correlation length is ^(t) = ^^r-^T^rliTfr + Rt) for 0<t<tq (or 
Tq <T <Tff.) and ^(t) = ^{Tg^p) when t>tq.liis assumed that at t = tq the temperature has 

a discontinuity while the domain size is continuous. 

If t* >tq, during cooling, the domains are in equilibrium. The domain size is 

r(0 = <§(r(0) for < r < , 

r(0 = min {r^ + a (Tg^p/T^) ln[l + (t- tq)/T(rq,Tg^p)\ , ^(Tg^p)] for t>tq , (7) 
with = .^(7^). 
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If <t* <tq, the domains are in equilibrium up to t* and fall out of equilibrium when 
t>t. Afterwards, they non isothermally grow from r(t*) = ^{Ttr + Rt*) down to Tq and 
then start to isothermally grow at T^^p. Therefore, the domain size is 



for 0<t<t 



r(t) = r(t*) + Uf ^^^^ '-^y for r < r < r„ , ( 8) 

* dt' 

t 

r(t) = min {r^ + a (Tg^p/Tfj) ln[l + (t- tq)/t(rq,Texp)] ' ^(Texp)} for ^>^q ' 



with rq = ^(Tq) 



For T„>T , the amplitude of aging of the domain size is 



' X 1 ^ 



For T >T„ and if T is not far from r„ , the same expansion as above reads 



^(Texp)-P(tq) = ^tr-^ 



f \ 
1 1 



y^^exp ^jTq J 



+ 



(^(Tq)-p(tq)) , 



with 



^(Tq)-p(tq) = ^t, 



T -T„ 



(ln(i?o/|i?|)f^^[5 + 31n(i?o/|i?|)] • (9) 



y a ) 



Consequently, for a given cooling rate R , there is a regime change at Tq = T . 

We emphasize that for the preceding two calculations the natural energy scale is not the 
microscopic barrier height T^y but the effective, much larger, barrier height f^. = T^j (^^^/a). 

Temperature cycling experiments are more complex. The comparison of the experiments 
with the predictions is only performed using numerical computation (see below). 
3.2. Deeper Insight into the Model - The previous simplified form of the model is able 
to describe some features actually observed in our experiments. However, it is not fully 
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satisfying for, at least, three reasons. First, the pictures of figure 5 clearly show that there is a 
distribution of domain sizes. This is not surprising if we know that the elementary growing 
process of a domain is a dipole flip: if a domain grows, one of its neighbors slims. Second, the 
effective barrier height T(j(ri) = (r^/a) Tg- may be seen as a mean value over all the frontier of 
a domain which crosses different regions of the sample where the density of pinning points is at 
random. Therefore, there must be a distribution of the height f(j . A third important point is that 
without a distribution of barrier heights, it is impossible to reproduce both the isothermal aging 
and the cooling rate dependence: either the (unique) barrier height is too large, and then there is 
cooling rate dependence but no aging at the final temperature, or the barrier height is too small 
and there is barely any cooling rate effects. Therefore, cooling rate effects must be attributed to 
large local barrier heights, which coexist with smaller barriers. Finally, our recent results on the 
response to a positive temperature jump from Ti to T2 (see sub-section 2-5) give again a very 
strong indication in favour of the coexistence of domains with very different sizes. Indeed, the 
non monotonic variation of e' cannot be found, even with a distribution of the domain sizes, if 
they all grow. The solution is to assume, immediately after the positive jump, that some 
domains grow because they are still smaller than ^(T2) coexisting with slimming domains 
because they have grown beyond ^(T2). 

The last point deserves some examination. We have seen that the delayed response of 
e' to a positive temperature jump is non monotonic and composed of a fast increase followed 
by a slow decrease. One might think that the increase of e' would be due to the decrease of 
domains by dipole flips at the frontier. This would be the inverse process of the growing 
process when the domains are large and consequently when the process is logarithmically slow. 
It would provide a very slow slimming process, in contradiction with the experiments. On the 
contrary, a fast process is provided by the birth (or nucleation) of new very small domains or 
nuclei inside the too large domains. These nuclei rapidly grow and finally induce the 
fragmentation of a large domain into new domains with their size equal to the coherence length 
^(T2). We first examine the fragmentation of a domain with a given size and afterwards we 
take into account the size distribution. 
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3.2.1. Fragmentation of a single domain - The process begins with the creation of a 
germ or nucleus which may appear anywhere in the too large domain. Its initial size is of the 
order of the lattice parameter a. It grows as 

p„(0 = a + a (T2/T„) ln[l + t/z(a,T2)] , 

but later its real size is 

r„(0 = min{a + a(r2/r^)ln[l + r/T(fl,r2)],.^(r2)} . (10) 

The index n stands for nucleus. For the present calculation we take as time origin t = the 
instant of the temperature jump. Some orders of magnitude have to be recalled. The size r(0) of 
a too large domain is such that ^{T2) ^ riO) <t,(Ti). Our experiments are performed with 
jumps from = 5 K to 72 < 20 K and consequently \ < ^{Ti)l ^{T2)<2. Therefore, a small 
number of nuclei growing up to ^{T2) are sufficient to break down the large domain into new 
ones with r„(0 < ^(T2). 

Let be r{T2) the number of nuclei created per volume unit and time unit at the 
temperature T2. If we take into account that new nuclei are continuously created, the 
contribution to the wall area built inside a given large domain, is, at short times, 

t 

(7^(0»r(r2)r(0)3 jdr'(r„(n)^ , (11) 


where the index A'^ stands for the whole set of nuclei inside the large domain of size r(0). We 
note that there is no factor of 1 /2 here because, at least at the beginning of the process, the 
nuclei are just like islands in a sea, without common coast. This equation is based on two 
assumptions: i) the size r(0) of the large domain is practically invariant during the 
fragmentation process; ii) the nuclei are created anywhere in the large domain at any time. The 
second assumption is unrealistic since new nuclei can only be bom in the part of the large 
volume not occupied by previously created nuclei. In fact, the nucleation process becomes less 
efficient as long as the number of nuclei increases. 
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At the beginning of the growing process t « T(a, T2) , when the second assumption is 
still satisfied, then p„(0 ~ a [l + (7'2/^<7) t/T(a,T2)\ and the newly created area in a single 

large domain is in lowest order 



aN(t) = r(T2)r(0f a^t 



(12) 



At ? = the new area and its time derivative are 
(7^(0) = and 



V dt Jt=Q 



^1 =r{T2)r{Qfa^ 



(13) 



The derivative depends on T2 only through r(T2). This is a fast initial increase of the area as 
we need. 

The process stops when the whole volume of the large domain is occupied by grown up 
nuclei. The condition which defines the time 0^ would read 



r{T2)r{0f jd?'(r„(n) «K0) 



(14) 



if the assumption ii) were true. Indeed, the value obtained in this way for 0^ is only a lower 
limit. From the definition of 0^, we have cr^(o°)~ <J^{6f^). Since o^{0) = 0, the total 
variation is simply Ao^ = Gp^(dj^) . 

Several cases have to be considered according to the relative efficiencies of the two 
processes (nucleation and growth) in action. This may be inferred from the first step of an 
iteration process which takes into account the remaining volume where new nuclei are allowed 
to appear. It gives 



(7N(t)'-r(T2)\df 



r(0)3-r(r2)r(0)3 jdr(r„(n)- 



For short times, this equation becomes 



<7Ar(0-a^r(r2)r(0)3 



t + 



^(T2/TW)__BM^^ 

T(a,T2) 2 



(15) 
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The initial straight line bends downwards if r(T2) a T(a,T2) (T(j/2T2)> I (strong 
nucleation) or bends upwards if r{T2) a T(a, T2) (T(j/2 T2) < 1 (fast growth). The physical 
meaning of the limit condition 



r{T2)a'' T{a,T2){T^llT2) = l 



(16) 



is that the number of nuclei created in a volume of size a during the characteristic time T(a, T2) 
of growth is of the order of unity. Indeed three cases have to be distinguished, separated by the 
two following limits. 

Nucleation is said to be very frequent if a new nucleus appears (in the mean) when the 
previous nucleus has hardly changed its size. The time 5t needed for doubling the nucleus size 
is defined by 



l + .^2 



5t 



= 2 a or 5t = T{a,T2) 
T2 



Tfj T(a,r2)y 

During this time, at least one nucleus must be bom. The limit condition reads 

r(T2)a^T(a,T2)(Ta/T2)^i ■ (17) 
It is calculated at the end of the fragmentation process, when the remaining volume is of the 
order of a . Within a factor of 2 the two equations (16) and (17) are identical. 

Nucleation is said to be very rare if a new nucleus appears (in the mean) only when the 
previous nucleus has reached its final size ^(T2). The time At needed for this evolution is 
defined by 



a In 



At 



T(j T(a,T2) 



= ^(T2) or At= z(a, T2) exp 



^(7^2) Tc 



'2 J 



During this time, at most one nucleus must be bom. The limit condition reads 



r(T2)r(0yt(a,T2)^W 



a T2 J 



<1 



(18) 



It is calculated at the beginning of the fragmentation process, when the remaining volume is of 
the order of r(0)^. 

The two previous limits define three nucleation regimes: 
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Regime 1 : rare nucleation 



r(T2)T(a,T2)< 



1 



r(0fcxTp{(^(T2)/a)(T^/T2)) ' 



Regime 2: intermediate nucleation 



1 



<r(T2)T(a,T2)< 



1 



r(0)\xp{(aT2)/a)(T^/T2)) 



a\T^/T2) 



Regime 3: frequent nucleation 



1 



<r(T2)T(a,T2) . 



a\T^IT2) 



The three regimes are examinated in detail in Appendix A. 

3.2.2. Distribution of the domain size - In an experiment with a positive temperature 
jump to T2 after a sojourn of duration ti at T^, the volumic fractions of large domains 

{^{T2)<r{0)<^{Ty)) and of small domains (r(0) < <^(r2)) are and = !-/'/«, 
respectively. The fraction Piaih ,T2) is an increasing function of T2 since when the 
temperature T2 increases the coherence length ^{T2) decreases and consequently, the interval 

^{T2)< r(0) <^{Ti) becomes broader. It is also an increasing function of which saturates 
for very long ti when nearly all the sample is occupied by domains larger than ^(T2). 
Moreover, its dependence on T2 is weakened when ti is long because most of the domains 
have size close to ^(Ti) while only a few of them are still in the size range of ^(72). 

Now we have to take into account the size distribution. Two cases must be 
distinguished: a) a simple experiment (cooling followed by isothermal aging) where all the 
domains play the same role; b) a more complex experiment (cooling, isothermal aging at Tj, 
positive temperature jump at T2 and finally isothermal aging at this temperature) where two 
types of domains must be separated: those larger than ^(T2) and those smaller than ^(T2). In 

the first case, any mean value is calculated through a sum over all the domains with the 
convenient weighting function. It is noted by brackets {F)^q^ with the index dom. In the 

second case, we have to add the contribution of the newly created walls through nucleation, and 
the final average reads: {F) - (F)^^ + Pia{F)]<^ where {F)j^ concerns only the walls of the 

nuclei. 
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Now we compare the results predicted by the model with our data. The second 
contribution in the above formula is well fitted by an exponential return 
^^incr [l ~ exp(-r/ Tj„^^)] (see sub-section 2.5) which has to be compared with the above term 
Pla {P)n which is proportional to P;^ (c^at);^. 

We first compare our experimental data on the amplitude ^Bincr t^e increase of e 
with the prediction /le^ oc p^^ (Ao^^)^^ provided by our model. The dependence on ti is 

obviously satisfied: the amplitude is an increasing function of ti which tends towards an 
asymptote when ti becomes very long. Concerning the dependence on T2, things are less 
clear. On the one hand, the observed amplitude ^e^^^^ decreases when T2 increases (see 
figure 3). On the other hand, in the expression for Ae the factor P/^ increases with T2 . 
Therefore, in order to have a decreasing product the factor I^Aoj^)^^^ must be decreasing too. 

This constraint selects the intermediate case and it additionally shows that the increase of the 
first factor is slower than the decrease of the second one. In this case, the amplitude Ae^^ is 



t 

(la[T2\n{d2lT^)-T^] 



^% - Pla ( .„r^ '-^"i'fwTJ ) ^^'^ 

la 



where 62 (defined as the time when the volume of a large domain is completely filled with 
newly created domains) is now given by the mean value 

^7- ^ ^ . (20) 

^{T2)a^[{T2|T^)\n{e2|T^)-\Y ) 

\ / la 

Provided that the dependence of r(T2) on temperature is less than Arrhenius, one indeed fmds 
that the second factor decreases as the temperature increases. Finally, we notice that the 
observed weakening of the T2 dependence when ti becomes very long is also predicted (from 
the behavior of P/^). 

Our second comparison is about the characteristic time. If we do the approximation 
the factor P;^ cancels and we get 
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1 



(21) 



ir{T2)a^ {[(T2lT^)\n{d2l 



o-i]) 



la 



Clearly, the two characteristic times T2 and G2 both decrease when T2 increases and they do 
not depend on ti. We recall that they are not equal because they correspond to different 
situations: the former is relative to the initial slope while the latter concerns the end of the 
nucleation process. 

Our experimental data are well fitted with the exponential law 
^^incr [1 ~ e^xp(-t/ Twiner)] • This means that the initial variation is governed by the time Tj^^^ 

which is found to be nearly independent of ti and to decrease when the temperature T2 
increases (see figure 4). We notice that the comparison with the characteristic times Ti, T2 
T3 of the three regimes (defined in Appendix A) is not a criterion of choice. Indeed, it does not 
rule out any regime since the three times vary like Tj„^^ (however, with the improbable 
additional condition for regime 3 that r(T2) 7^2 '^^^ decreases when temperature increases). 

As a conclusion, we emphasize that the main features of the non monotonic behavior 
which immediately follows a positive temperature jump are quite well accounted for by our 
modified (and improved) model. 

3.3. Numerical Calculations on the Model - The different behaviors of the domain 

size r(t) predicted by our model in its simplified form were computed for a unique barrier 
height and consequently no size dispersion. We recall that in this case the variation of the 
dielectric constant is simply proportional to the inverse of the unique domain size K /r(t) and 
its lowest limit \% K i^exp ) • 

The results concerning r{t) are sketched in figure 6 which shows the reduced domain 
size r(t)l^(Tg^p) (full line) for the barrier height and a cooling from temperature Tf^ to 

temperature T^^p lasting the time = (r^^ - T'g_^p)/|^| followed by an isothermal evolution. 
The data are displayed together with a reference. The ascending curve 



^{m)/^(T,^) = (^{Ttr)/^(T,^)) ^T^/(Ttr + Rt) shows the upper limit for r(t)/^(T,^p) 



(dashed line). If T^j is large, the size cannot increase and r(t) = at any time. If T^- is very 
small, the domains are always in equilibrium and r{t) = ^{T{t)) until t = t(. and then 
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r(t) = ^{Tg-^p). For intemediate values of T^j such as that used in figure 6, the representative 
point also begins with a travel on the curve corresponding to r{t) = ^{T{t)) until t = t* (or 
T = T*) where it leaves the curve, then goes slowly upwards and after t = t(., logarithmically 

increases. Two different time scales are used in the two half-figures in order to emphasize these 
features: a) the left half-figure covers the cooling time t^. and the splitting between r(t) and 
^{T(t)) is clearly seen at time t* or temperature T* ; b) the right half-figure covers 10^ x and 
the very slow isothermal growth of r{t)l ^{T^-^p) becomes obvious. This illustrates very well 

the set of equations (3). 

Indeed, from the experimentalist's point of view, it is more convenient to take as time 

origin the instant when isothermal aging begins. Hence, the relevant time is either t = t-t(. 
(full ramp) or t = t-tq (quenching). In the experiments, the variable t is indeed t . In what 

follows we adopt this new time origin. Here we present the results of our computations on the 

time dependent part of the dielectric constant (in arbitrary units) during the isothermal evolution 
following a temperature ramp in the form 5e{t) = K {ijrit) -Ij^iTg^p)^. In this way, if the 

domains have reached their equilibrium size ^(Tg^p), then 5e = 0. Our data are displayed in 

Table I. We have calculated the values of d£(i) at f = 0s (beginning of the isothermal 
evolution) and f = 10^ s (a duration much longer than our experiments) for different values of 
the unique effective barrier height r^- = T^- (^ff/a). We have used for the cooling rate R two 

extreme values which span the range of our experimental values. 



Table I 



These results deserve some comments. For T^j = 40 K the domains always have their 
equilibrium size: there is no aging. For Tq- = 60 K the domains have not their equilibrium size 

Q 

at the end of the ramp, but they are able to reach it in less than 10 s: this is ergodic aging. For 
r^j- = 90 K and r^j- = 100 K the domains have not their equilibrium size at the end of the ramp 
and they are not able to reach it within 10 s: this is aging with effective ergodicity breaking 
(effective, since the model is fundamentally ergodic). For T(j = 150 K the domains have not 
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their equilibrium size at the end of the ramp and the growing process is so slow that it cannot be 

observed during 10 s: this looks like the case of no aging, except for the dependence on R. 

Two other features can be seen in Table I. Firstly, the maximum possible evolution 
^^max ~ ^^(0) - 5e(oo) = 5e(0) monotonously increases with . Secondly, the evolution 

Q 

observed during a long time Sei^^g = d£(0) - d£(lO s) is equal to for both low barriers and 
high barriers, and passes through a maximum for in the middle range. This is easy to 
explain. For small the domains rapidly grow, they are able to reach the coherence length 
within a reasonable time, dsiong equal to 5e^^, and they both tend to if T^j is very small. 

For large the domains can hardly grow, during the ramp they shortly leave the coherence 
length, they cannot reach the equilibrium in a finite time, and ^C/^^^ is smaller than Se^^. If 
f(j is very large, they do not grow at all at T^^p, and therefore 5£/on^ =0 - For intermediate 
T(j values, 5e/^„^ which starts from and goes to 0, passes through a maximum. 

These computations provide other interesting informations. We call e^j„ the value 
obtained after 400000 s. This is the duration of our longest experiments. Then we have 
shown: i) e^jyj is nearly (but not perfectly) a linear function of the logarithm of the cooling rate, 
at least for 0.001 K.s~^ < < 1 K.s~^; ii) the cooling rate R and the relaxation time exactly 
play the same role (indeed, e^;„ is a function of the product R only). 

Numerical calculations were also performed when there is a sudden change from 
temperature to temperature T^-^. The dependence of e^j„ on the quenching temperature 

is strong for close to T^^. and weak at low temperatures. Temperature cycles were also 
simulated. 

Finally, we have seen that the simplified form of our model provides a rather good 
qualitative agreement with our data for most of our experiments: simple aging, quenching, 
positive or negative temperature cycles. However, as discussed above, some qualitative 
features (such as the coexistence of cooling rate effects and isothermal aging) require some 
distribution of the domain size. We try now to deduce more information on this distribution, 
and find interesting qualitative conclusions. 
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4. Distribution of barrier heights 



In a series of experiments all performed with the same cooling rate i? = -0.015K.s ^ 
on a KLT sample with the lithium concentration x = 0.017, we have observed that the 
amplitude of the decay de= ^£(0) - <5e(10^ s) during 10000 s varies in a non monotonic 
manner with the temperature T^^p of the isothermal aging. The data are displayed as full 

squares in figure 7. They present a broad maximum near 7 = 15 K. It is actually impossible to 

reproduce this feature with a single barrier height model. Rather, we attempt to determine the 
full distribution function P(T(j) of the barrier heights using this effect. We choose to 

parametrize it using a gamma distribution which contains only two free parameters and gives a 
rather good fit. This normalized function reads 

Y(T ) fT f T ^ 

P(Ta)= — exp with v>0 , (23) 



r(v)0 









where Y(T(j) is the Heaviside step function and r(v) is the gamma function. The distribution 
is scaled by the parameter 0, related to the most probable barrier by = (v-l)0if V>1; 
if 0< v<lthen rj^ = 0. 

In figure 7 are shown the data together with the results of the calculation with three 
values of the exponent v. For v = 1 (then P(T(j) is simply a decreasing exponential) the curve 

steadily decreases, in clear contradiction with the experiments. The curves for the parameter 
= 60 K and the exponents V = 2.2 and V = 3 both present a maximum in the right place. 
However, for v = 3 the curvature is too sharp, while for v = 2.2 the fit to the data is quite 

good. Moreover, we have established that the position of the maximum of the curve roughly 
varies as the most probable value T^^ . These results mean that the data cannot be reproduced 

by a distribution which does not vanish for small barriers. More precisely, a good fit needs that 
V = 2, which implies that the distribution varies as P{T(y)^T(y for T(j«0, and that 
r^^ = 60K. We now adopt these values. Then the corresponding mean value is 
jmean ^ ^qq ^^^^ ^-^j^j^ ^^^^ hQighi is AT^j = 155 K. In principle, the result is 

valid only for the x = 0.017 sample, but we assume that it retains the characteristics of all our 
KLT samples. We therefore use the distribution 

26 




(23) 



The other needed numerical quantities are the lattice parameter a = 4 x 10 m, the coherence 

length at the transition temperature ^(Tf^) = 5 a and the relaxation time T^o - 10~^^ s. The 
computed quantity is simply the ratio of mean values [p^'^ j {p^^ the variation of which is 

proportional to 5e{t) , normalized at the largest value. 

We only report our numerical calculations for the very important results concerning the 
dependence of the asymptotic value of e(t) as the limit of the decay following a temperature 
ramp at the constant cooling rate R. The experiments clearly show a logarithmic dependence 
[8], also predicted by the model. The results from the model are displayed in figure 8. They 
show that the asymptotic value depends on the cooling rate as e = p + q log\R\ at least in the 

range from - 0.001 K.s~^ to - 0.1 K.s~\ as in our experiments. The logarithmic variation is 

perfectly obeyed, actually more accurately than it is when the curve is calculated without 

distribution of the barrier height. This is again very satisfactory. 

Moreover, we have examined the role of the barrier distribution on the coefficient q of 
the logarithm. We have found that q depends on the unique parameter of the distribution 

in a non monotonic manner: the coefficient q{T^^) passes through a maximum for 

T^^ = 110 K. Our experiments [8] show that q increases with the lithium concentration x 

when it is smaller than = 0.02 and the Li+ ions freeze at into the dipolar glass phase. 

Indeed, when x increases, the density of pinning points in the crystal increases too, leading to 
larger equivalent barriers and consequently to a larger T'^^. This is in agreement with our 

computations. When the concentration x is larger than x^ =0.02 the Li+ ions undergo a 
transition at Tf^. towards the ferroelectric phase and we observe that the coefficient q then takes 
small values [8]. We attribute this decrease not to a low T'^^ value but on the contrary to a 
large one, larger than = 110 K. This is probably due to the correlations between the 

random fields which must exist in this phase. Then the large density of correlated pinning 
points is equivalent to a barrier distribution shifted towards high energies. 

5. Conclusion 
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In this paper we present an extended version of a brief report [10] on the model 
elaborated in order to explain a large series of experiments performed on a family of disordered 
dielectric crystals. The central idea of the model is to attribute the slow isothermal decay of the 
alternative dielectric constant e(co,t) to the motion of the walls of the polarized regions or 
domains of the crystal lattice. This motion is extremely slow because it is hindered by the static 
random fields associated with the lithium ions which are frozen at the experiment temperatures; 
moreover it becomes slower and slower as long as its size increases. Thus the existence of such 
long-lived polarization regions in the paraelectric phase allows us to use the words domain and 
wall. 

In addition to well known features of slow dynamics and aging, such as effective 
ergodicity breaking with logarithmic dependence on the cooling rate, we report new 
experiments on a rather surprising non monotonic response of e to positive temperature jumps. 

Along the lines of the model we calculate the evolution of the typical domain size, 
assuming in a first step that the growth is slowed down by the same barrier height for all the 
domains. Even with this oversimplified assumption we are able to qualitatively account for 
many of the observed effects. 

However, the very coexistence of cooling rate effects and isothermal aging, the non 
monotonic time dependent response to positive temperature jumps and the non monotonic 
temperature dependence of the amplitude of isothermal decay can be explained only if a broad 
distribution of the domain barrier height is taken into account. It is satisfying that these three 
apparently uncorrected features lead to the same conclusion. Therefore, in a second step we 
have computed some typical behaviors taking into account the distribution of barrier height T^j ; 
we have shown that it goes to with 7^. and, consequently, that it has a maximum for a non 
zero value of T^y . 

As a matter of fact, if a broad barrier distribution is necessary to explain the non 
monotonic response to a temperature jump, it is not sufficient to induce the observed initial 
rapid variation. For that purpose, we have introduced a nucleation process which is known to 
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be fast because it concerns small entities. The results are then in good qualitative agreement 
with experiments. 

Finally, we emphasize that the model is not very complex since the dynamics of the 
domains is completely characterized by a single time dependent length r and the function which 
describes the barrier distribution, itself parametrized by its most probable value. Nevertheless, 
this model is quite powerful since it describes semi-quantitatively all the features observed in 
KLT crystals. 
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APPENDIX A 



The three different nucleation regimes are now considered, taking into account the decay 
of the remaining volume in the too large domain. 

Regime 1 - If nucleation is very strong while growth is slow, the process stops at 
- when the total volume of the large domain is filled by very small domains of typical 
size a which had not enough time to grow. Under this oversimplified hypothesis, the 
unrealistic assumption that the nuclei are created anywhere in the large domain at any time can 
be discarded. Then the evolution of the number N(t) of nuclei is governed by the equation 



— = r(T2)a 
at 



rioy 



-N 



V 



the solution of which is 



N(t) = ^ (l - exp(-r(r2) t)) 



Consequently, the total nucleus area is 
r(0)3 



c^Ar(0 = — (l-exp(-r(r2)a^)) . 



The expansion of equation (A-1) for short times leads to 



(7^(0 = r(r2)r(0)3 



V 



(A-1) 



As expected, it is identical to equation (15) with x{a,T'£) = oo. With = l/(r(T2) a^), the 
fragmentation process is practically achieved for 0^ » 10 . 

Indeed, the final total area is c7jy(oo) ~ r(0)^/(2 a) since a factor of 1/2 has to be 

included since all the domains are in close contact and consequently the total variation is 
Ao^ ~ r(0)^/(2 a). Therefore, we adopt as a useful approximation the equation 



(^Nit) = — fit) (l - exp(- r(T2) t)) 



(A-2) 
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The slowly decreasing function (p(t) varies from 1 for f = (isolated nuclei) to 1/2 for ? — > oo 
(nuclei in close contact). We notice that AOf^ does not depend on T2 . 

In this case, the new domains have still to grow towards the equilibrium size ^(T2) for 
t> Oi. This is the standard aging process which induces a decrease of the dielectric constant. 

Regime 3 - If growth is very fast and nucleation is very slow, when a nucleus was bom 
attime?',its size is r„(0 = min {a + a (72/7^^) ln[l + (?- 0/'Z^(a>72)] > "^(^2)} at time t>f 
and it reaches its final size ^(72) after the delay 
At = t-t' ~ T(a,T2) Q^p\{^(T2)/a) [T(j /T2)\. The process stops at 0^ = ^3 when the total 

volume of the large domain is filled by domains of typical size ^(T2). Then the evolution is 
exactly the same as in regime 1 with a changed into ^(72). Here the average number of such 
new domains in their host domain is equal to (r(0)/^(T2)) . Consequently, the final total area 
is t7;^(oo) = r(0)Y(2(^(r2)). The total variation Aoj^ = r(Of /{2 ^(T2)) is an increasing 
function of T2 and the area evolution is governed by the equation 

(JN(t) = ^^m(i-^W(-nT2)^(T2ft)) . (A-3) 
S(^2) 

With T3 = l/(r(T2) ^(T2)^), the fragmentation process is practically achieved for % = 10 T3. 
In this case the new domains have reached the equilibrium size ^(T2) and they no longer 
evolve after ^3 . 

Regime 2 - In the two previous regimes when a new nucleus appears the already 
present nuclei are not evolving: in regime 1 because growth is extremely slow; in regime 3 
because growth is achieved. On the contrary, in the intermediate case when a new nucleus was 
bom the present nuclei are still growing; indeed most of them have not yet attained the 
equilibrium size ^(T2) and they grow further. 

If the assumption that the nuclei are created anywhere in the large domain at any time is 
maintained the characteristic time 0^ = ^2 is determined by the equation 

O2 

r(T2)a^ \dt' ((T2/T^)ln[f /T(a,T2)]f '-I . 
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As already pointed out, the logarithm is a slowly varying function. Consequently, this function 

may be taken as constant in the integral, with t' ~ 62- This approximation, together with 
t(a,T2) = Too exp( 7(^/72), leads to 

r(T2)a^ d2[(T2/T^)\n(d2/T^)-lf . (A-4) 

As usual, the value of 62 in the logarithm is not important and equation (A-4) provides a lower 
limit for 62 which reads 

62- ^y-^ 7j . (A-5) 

r(T2)a^ [T2ln(d2 /t^)-T^Y 

The creation of nuclei is more efficient if T2 is higher since the rate r(T2) is an increasing 
function of T2. Moreover, it is independent of the enclosing domain size r(0). Therefore, 
equation (A-5) shows that 62 becomes shorter when T2 is increased and that it depends 
neither on the time ti elapsed at Ti nor on the size r(0) of the large domain. In the present case 

the average number of such new domains in their host domain is equal to 

[r(0)/((r2/r(j) a In[02/T(a,r2)])j . Consequently, the final total area is 

a^(oo) = r(0)^/(2 {T2/T(j) a In [02 /T(a, T2)]) . The total variation is 

AaN = 7 — r , (A-6) 

2a{(T2/T^)ln(d2/too)-l) 

a decreasing function of T2. The stopping condition implies that p„(02) < ^(7'2)- Therefore, 

when t> 02 the new domains still have to grow towards the equilibrium size ^(T2). This is 

the standard aging process which induces the decrease of the dielectric constant. 

The time variation is expected to be of the form (7^(0 = Aop^ f{t) with /(O) = 0, 

/(oo) = 1 and f{t) oc t for small t. A convenient function is f(t) = 1 - exp (- ? / T2 ) • Then 



V dt Jt=Q T2 
and with equations (A-5) and (A-6) it comes 
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According to this equation, the characteristic time T2 does not depend on the time ti elapsed at 
Ti but it may depend on T2 . 

Our experiments show that the total increase ^ctjv is smaller when T2 is higher. This 
feature necessarily eliminates the first two cases but keeps the last one (see sub-section 3.2.2). 
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TABLE I 





R = 


- 0.001 K.s~^ 




-0.1K.s~^ 




i = Os 


f = 10^s 


f = 0s 


f = 10^s 


40 














60 


6 





17 





90 


33 


26 


45 


28 


100 


41 


40 


54 


49 


150 


73 


73 


88 


88 



TABLE CAPTION 

Table I - Computation of the time dependent part d£(t) = K J r(t) - 1 J ^(Tg.^)j of the 

dielectric constant at the beginning ( f = s) of the isothennal evolution following a temperature 
ramp and after a long time (t =10 s). A single domain size r(t) is assumed, the growth of 
which is characterized by the unique effective barrier height f^- = T^- (^^^/a). Two extreme 

values of the cooling rate R are used. The data are given in arbitrary units. 
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FIGURE CAPTIONS 



Figure 1 - Variations of the capacitance C after positive temperature jumps from 7^ = 4.3 K 
to several higher temperatures T2. From top to bottom: 11.3 K; 72 = 28.9 K; 

72 = 33.8 K; 72 = 39.2 K. The transition temperature of the sample is 7j^=35K, 
corresponding to x = 0.017. The measurement frequency was / = 100 kHz. Before the jump, 
the sojourn at 7^ lasted ti = 5000 s. For clarity the curves are vertically shifted relative to each 
other. 

Figure 2 - Variations of the capacitance C after positive temperature jumps from 7^ = 4.3 K 
to 72 = 16.2 K after several durations ti of the sojourn at 7^. From top to bottom: ti = 500 s, 
ti = 10000 s and ti = 80000 s. The transition temperature of the sample is 7^^ = 35 K. The 
measurement frequency was / = 100 kHz. The temperature jump times are taken as the 
common time origin. 

Figure 3 - Amplitude ASincr of the variation Ae^fi^^ [l- exp(-?/Tj„c^)] of the dielectric 
constant after positive temperature jumps as a function of the duration ti of the sojourn at 7^. 
The temperature jumps are from 7^ = 4.3 K to 72 = 16.2 K (squares), 72 = 18.7 K (circles) 
and 72 = 23.8K (triangles). The transition temperature of the sample is Tfj. =35K. The 
measurement frequency was / = 100 kHz. The data are displayed on logarithmic scales. The 
amplitude is given in arbitrary units. 

Figure 4 - Characteristic time Tj>j(,^ of the variation /lej„^^ [l-exp(-r/Tj>j(,^)] of the 
dielectric constant after positive temperature jumps as a function of the duration ti of the 
sojourn at T^. The temperature jumps are from 7i=4.3K to 72= 16.2 K (squares), 
72 = 18.7 K (circles) and 72 = 23.8 K (triangles). The transition temperature of the sample is 
Tf^ = 35 K. The measurement frequency was / = 100 kHz. 
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Figure 5 - Monte Carlo simulation of the three-dimensional random field Ising model at 
constant temperature T <Tg. The figures are snapshots of a cross-section of the system, taken 

at four times < ^2 < ^3 < U after preparation at t = 0. The initial condition is random Si=±l 

with probability 1/2 and corresponds to an infinitely rapid quench from T = oo to T which has 

been performed at t = 0. 

Figure 6 - Reduced domain size r(ty^(Tg^p) (full line) for a unique barrier height 
= 18 K and a cooling from Tf^ = 35 K to Tg^p = 5 K at the rate = - 0.01 K.s~^ followed 
by an isothermal evolution. The reduced coherence length ^{T(t))l^(Tg^p) is also shown 

(dashed Une). Using two different time scales, two important features are put in evidence: a) the 
left half- figure covers 3000 s and the time t* where r(t) and ^(7^(0) split off from each other 

is clearly seen; b) the right half-figure covers 3x10 s and the very slow isothermal growth of 
r(ty^(Tg^p) becomes obvious. 

Figure 7 - Amplitude of the isothermal decay of the capacitance C (normalized by the largest 
value) during 10000 s, measured at / = 100 kHz, as a function of the temperature Tg^p after 

the same cooling rate R = - 0.015 K.s~^ The data points (full squares) were obtained on the 
sample with x = 0.017. The three curves are calculated with the gamma distribution (see text). 
The exponents are v = 1.0 (dotted line), v = 2.2 (full line) and v = 3 (dashed line), 
respectively. The best fit is obtained for v = 2.2. 

Figure 8 - The dielectric constant (in arbitrary units) calculated with the gamma distribution 

(exponent V = 2, most probable value at 60 K) after an isothermal decay lasting 400000 s at 
r = 5 K and displayed as a function of the cooling rate on a semi logarithmic scale. The 
calculated points exactly fall on a straight Une. 



36 



REFERENCES 



[I] PALMER R. G., Adv. Phys. 31 (1982) 669. 

[2] BOUCHAUD J.-P., /. Phys. I France 2 (1992) 1705. 
[3] HOCHLI U. T., KNORR K. and LOIDL A., Adv. Phys. 39 (1990) 405. 
[4] VUGMEISTER B. E. and GLINCHUK M., Rev. Mod Phys. 62 (1990) 993. 
[5] KLEEMANN W., SCHONKNECHT V. SOMMER D. and RYTZ D., Phys. Rev. 
Lett. 66 (1991) 762. 

[6] DOUSSINEAU P., LEVELUT A. and ZIOLKIEWICZ S., Europhys. Lett. 33 (1996) 
539. 

[7] DOUSSINEAU P., LEVELUT A. and ZIOLKIEWICZ S., Europhys. Lett. 33 (1996) 
391. 

[8] ALBERICI P., DOUSSINEAU P., and LEVELUT A., J. Phys. I France 7 (1997) 
329. 

[9] ALBERICI-KIOUS P., DOUSSINEAU P., and LEVELUT A., Europhys. Lett. 39 
(1997) 329. 

[10] ALBERICI-KIOUS P., BOUCHAUD J.-P., CUGLIANDOLO L. P., DOUSSINEAU 
P. and LEVELUT A., Phys. Rev. Lett. 81 (1998) 4987. 

[II] COMES R. and SHIRANE G., Phys. Rev. B 5 (1972) 1886. 
[12] SCOTT J. P., Rev. Mod Phys. 46 (1974) 83. 

[13] VAN DER KLINK J. J., RYTZ D., BORSA P. and HOCHLI U. T., Phys. Rev. B 27 
(1983) 89. 

[14] VINCENT E., HAMMANN J., OCIO M., BOUCHAUD J. P., CUGLL\NDOLO L. 
P., in Sitges Conference on Glassy Materials, edited by M. Rubi (Springer- Verlag, 
Berlin, 1997). 

[15] HOCHLI U. T., WEIBEL H. E. and BOATNER L. A., J. Phys. C 12 (1979) L563. 
[16] STACCHIOTTI M. G. and MIGONI R. L., /. Phys. Cond. Matter 2 (1990) 4341. 
[17] NATTERMANN T., Ferroelectrics. 104(1990) 171. 



37 



figure 1 




38 



figure 2 




39 



[ure 3 



10 - 




100 



1000 



10000 



100000 



40 



figure 4 



2000 



1500 - 



1000 - 



500 - 




100 



1000 



10000 



100000 



41 



figure 5 



figure 5a 




42 



figure 5c 




43 



figure 6 




44 



[ure 7 




45 



[ure 8 




46 



